arXiv: 1503.0817lvl [math-ph] 27 Mar 2015 


On the integrability of a system describing the stationary 
solutions in Bose-Fermi mixtures 

Ognyan ChristoF and Georgi Georgiev^ 
t Faculty of Mathematics and Informatics, Sofia University, 

1164 Sofia, Bulgaria 

1 Department of Mathematics and Informatics, University of Transport, 

1574, Sofia, Bulgaria 


Abstract 
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1 Introduction 

In this paper we study the integrability of the system that comes from the time dependent 
mean held equations of Bose-Fermi mixture (BFM) in one dimensional optical lattices. The 
interest in BFM arises after the discovery of Bose-Einstein Condensates (BEC) in 1995 and the 
desire to understand strongly interacting and strongly correlated systems, with applications in 
solid state physics, nuclear physics, astrophysics, quantum computing and nanotechnologies. 
For more detailed physical background of BFM we refer to p in m ei si and the literature 
therein. 

At mean held approximation we consider the following Nf -(-1 coupled nonlinear Schrodinger 
equations 
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where the wavefunctions dU describe each of Nf fermions and d/ b is the wavefunction for the 

Nf 

bosonic component, pj = ]U |dF| 2 and <?bb, <?bf, ^f, itib are certain physical constants. In 

2—1 


( 1 . 1 ) 

( 1 . 2 ) 


1 




particular, g B B and g BB are related with the s-wave collisions for boson-boson and boson- 
fermion interactions, respectively. The potential V is usually of the form V = Vosn 2 (ax,n), 
where sn(ax, k) is the Jacobi elliptic sine function. In this paper we take Vo — 0 as in [1]. 
We are interested in the stationary solutions to the system (11.11) . (11.21) of the kind 
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where n 0 , kqj are constant phases, qo,qj and 0 O , Qj are real-valued functions related by 
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C Q , C 3 , being constants of integration. After substituting (ll.3|) . (jl.4[) in equations (II.ip . (II . 2[) 
and separating the real and imaginary part we get 
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Kostov et al. [9] have found plenty of particular (quasiperiodic, periodic and soliton) 
solutions to the system (11.61) and therefore, stationary solutions to the system ( 11 . 11 ) . ( 1 1 . 2 j) . It 
is natural to ask whether we can obtain more, that is, for what set of constants the system 
(II. 6 p has enough first integrals to be integrable. Note that when g BB = 0 the equations 
separate, i.e., the system is solvable. 

Before giving our main result let first get rid of the inessential (for integrability) param¬ 
eters. In what follows we assume that the parameters cuo, ujj, mp, me, g BB are positive since 
they have an origin from physics, and Cq, Cj, g B F are arbitrary real parameters. We put 
% = f3q 0 , qj = otijj, x = jx. Then we choose a = v /mF, /3 = y/rri^, 7 = l/(m By /g^), g BB ^ 0. 
Denoting g BF = 9BFO/ 2 ^ 2 m B , D 0 = ^o7 2 ? ™b,<A = uJj^ 2 m F ,C 2 = C 2 ^ 2 /a A ,Cl = C^ 2 / (3 A we 
reach 

1 d 2 q 0 _ 3 _ (_ 2 \ _ _ _ Cq 

2~dSC ~ % ~ 9B¥ l ^ Qi 1 Qo + UoQo ~ 2^’ ( 1 ‘ 7 ) 

1 d 2 qj „ _ 2 ~ , ~ ~ Cj , N 

9BF%Qj + UjQj = W J = 1} ''' ’ Nf ' 

To simplify notations we skip the tildas, write t instead of x and denote Pj — <jj, j — 0,..., Nf, 

(. = d/dt ). Then the system (1 1. 7[) can be presented as a Hamiltonian system with the 

Hamiltonian 
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For the Hamiltonian system with the Hamiltonian (II. 8 p we consider the cases: 

1) Co = 0 ,Cj + 0 : u)j = u 2 /2,j = 1,,7V/ ; 

2) C 0 ^ 0,Cj — 0,j — 1,..., N f ; 

3) Co 7 ^ 0, Ci 7 ^ 0, Nf = 1„ g-BF sufficiently small. 

Our result is the following: 

Theorem 1. For the cases given above, the Hamiltonian system corresponding to M.8\) is 
non-integrable in Liouville sense unless g-BF = 0. 

In other words, the Hamiltonian system under consideration is integrable only when it is 
separable. 

The proof of the above result is based on the Differential Galois approach and Ziglin- 
Morales-Ramis method. This method has been applied for the studying the integrability 
to a number of Hamiltonian systems, in particular systems with homogeneous potentials, 
see The classification of all integrable two degrees of freedom systems with 

polynomial potentials of degree 3 is obtained in mu- In particular, the above mentioned 
approach is used in [ 1 ] for obtaining non-integrability results for some two degrees of freedom 
Hamiltonians with rational potentials. Note that the system in this paper is not of that kind. 
For the natural Hamiltonian systems with two degrees of freedom, similar to (II. 8 p 

H = dM + U(guq2) 

there is an integrable generalization of Garnier’s system found by Wojciechowski [2lj, namely 

U = Aq\ + Bql + {q\ + q 2 2 f + % + 

Qi Q2 

with a rational first integral depending on A, B, C, D (see also [16]). Note that in the system 
under consideration, the symmetry is lost, so it is natural to expect integrability only in the 
separable case. 

The paper is organized as follows. In the next section we recall some facts about Differ¬ 
ential Galois groups and Morales-Ramis method which we use. Then, Section 3 is devoted to 
the proof of Theorem 1. We finish with some comments. 


2 Differential Galois Theory and Integrability 

Here we summarize some notions and results related to Ziglin-Morales-Ramis theory. 

A differential held is a held with derivation d — , i.e. an additive mapping satisfying 
Leibnitz rule. A differential automorphism of K is an automorphism commuting with the 
derivation. 

Consider a linear system 

x = A[t)x , x E C n (2.1) 

with t dehned on some Riemann surface. Denote the coefficient held in (12.ip by K. Let x t] be 
the elements of the fundamental matrix X(t). Let L(xij) be the extension of K generated by 
K and Xij - a differential held. This extension is called Picard-Vessiot extension. Similarly 
to classical Galois Theory we define the Galois group G := Galx^L) = Gal(L/K ) to be the 
group of all differential automorphisms of L leaving the elements of K hxed. The Galois 
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group is, in fact, an algebraic group. It has a unique connected component G° which contains 
the identity and which is a normal subgroup of finite index. The Galois group G can be 
represented as an algebraic linear subgroup of GL(n, C) by 

v(X(t)) = X(t)R a , 

where a G G and R a G GLfn , C) (see e.g. [2Dj) - 
Consider now a Hamiltonian system 

x = Xh(x), t G C, x G M (2.2) 

corresponding to an analytic Hamiltonian H, defined on the complex 2n-dimensional manifold 
M. Suppose the system (\2.2\i has a non-equilibrium solution T(f). Denote by T its phase 
curve. We can write the equation in variation (VE) along this solution 

£ = DX H {*{t))Z, f G T r M. (2.3) 

Further, using the integral dH we can reduce the variational equation. Consider the normal 
bundle of T, F := TyM/TM and let tt : TyM —» F be the natural projection. The equation 
(12.3p induces an equation on F 

f] = vr,( J DX H (T(t))(7r- 1 7 / ), rj G F. (2.4) 

which is called the normal variational equation (NVE). 

It is natural to assume that if the system (12.2p is integrable, then the linear equations 
(VE) and (NVE) are also integrable. 

The solutions of (12.3p define an extension Ly of the coefficient field K of (VE). This 
naturally defines a differential Galois group G = Gal(Ly/K). Then, the following result has 
established 

Theorem 2. (Morales-Ruiz-Ramis |12] ) Suppose that a Hamiltonian system has n mero- 
morphic first integrals in involution. Then the identity component G° of the Galois group 
G = Gal(Li/K) is abelian. 

Once it is proven, that G° is not abelian, the respective Hamiltonian system is non- 
integrable in the Liouvillc sense. Note that the fact that G° is abelian doesn’t imply nec¬ 
essarily integrability of the Hamiltonian system. Thus, one needs other obstructions to the 
integrability. A method based on the higher variational equations has been introduced in [12] 
and the previous Theorem has been extended in [T3] . Before formulating this result let us 
give an idea of higher variational equations. For the system (12. 2 p with a particular solution 
'F(f) we put 

x = T(f) + &£ {1) + e 2 f {2) + ... + e k ^ k) + ..., (2.5) 

where e is a formal small parameter. Substituting the above expression into Eq. (12.2p and com¬ 
paring terms with the same order in e we obtain the following chain of linear non-homogeneous 
equations 

£ {k) = A(t)^ + / fc (£ (1) , • • •, k = 1, 2 ,..., (2.6) 

where Aft) = DXni^ft)) and f\ = 0. The equation (12.6K is called k-th variational equation 
(VEfc). Let Xft) be the fundamental matrix of (VEi) 

X = Aff)X. 
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Then the solutions of (VE*,), k > 1 can be found by 

^ = X(t)c(t), (2.7) 

where c(t) is a solution of 

c = X~\t)f k . (2.8) 

Although (VEfc) are not actually homogeneous equations, they can be put in that frame, and 
therefore, one can define successive extensions K C L\ C L 2 C ... C L k , where L k is the 
extension obtained by adjoining the solutions of (VE*.). Correspondingly one can define the 
Galois groups Gal(Li/K ),..., Gal(L k /K). The following result is proven in [1 13], 

Theorem 3. If the Hamiltonian system Ii2.2\) is integrable in Liouville sense then the identity 
component of every Galois group Gal(L k /K ) is abelian. 

Note that we apply Theorem [3] in the situation when the identity component of the Galois 
group Gal(Li/K) is abelian. This means that the first variational equation is solvable. Once 
we have the solution of (VEi), then the solutions of (VE*,) can be found by the method of 
variations of constants as explained above. Hence, the Galois groups Gal(L k /K) are solvable. 
One possible way to show that some of them is not commutative is to find a logarithmic term 
in the corresponding solution (see detailed descriptions and explanations in [121 13l [T3]). 

Now we recall a perturbational technique which is still related to the Differential Galois 
approach. Let Mo be a two-dimensional complex analytic symplectic manifold, Ho(q,p) be a 
holomorphic Hamiltonian and Xh q be the corresponding Hamiltonian vector field. Assume 
that the system 

q = H 0tP , p = -H Q>q (2.9) 

has a hyperbolic equilibrium ( qo,Po )• Then the system (12.9p has a separatrix 

To : (qo(t),Po{t)), lim q 0 (t) = q 0 , lim p 0 {t) = p 0 - (2.10) 

t—> oo t—> oo 

The functions qo(t),po(t ) are meromorphic in f 6 C. Let 

H(q,p,t,e) = H 0 (q,p) +eH 1 (q,p,t) + ... ( 2 . 11 ) 

be a meromorphic small (complex) perturbation of H 0 satisfying Hi(q,p,t + oj) = Hi(q,p,t) 
with a period oj G C. This function H is defined over M = M 0 x F u , F u = C/cnZ. We can 
write the Hamiltonian system defined by H(q,p,<p) over M as 

q = H p , p=-H q , 0 = 1, (q,p,(p)eM. (2.12) 

When e = 0 the system (j2 .1 2j) reduces to 

g = M 0 , p , p = -H 0>q , 0 = 1, {q,p,ip)eM. (2.13) 

The unperturbed system (12.13jl has a hyperbolic cu-periodic orbit n 0 := (qo,po,ip = t( 
mod w)). It is well known that for small | e\ the perturbed system ( 12 . 121 ) has also an u- 
periodic orbit n e := (q(t, e),p(t, e),<p — t — t 0 ( mod u>)), such that (q(t, 0 ),p(t, 0 )) = ( qojPo )• 
We define the (stable) complex separatrix A0 of the system (12.121) as the set of integral 
curves of (12.121) asymptotic to n e as t —y oo. For fixed e, it is a two-dimensional complex 
surface. This separatrix can have transverse self-intersections. 
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Remark 1. Recall that in the real case the separatrices can not have transverse self¬ 
intersections. Such intersections can occur between stable and unstable separatrices. For 
real Hamiltonian systems, the existence of such transverse orbits is considered as a source of 
chaotic behavior and is an obstruction to existence of an analytic first integral. 

Ziglin [22] proved that for complex Hamiltonian systems, the existence of transverse self¬ 
intersections for separatrices is also an obstruction to the integrability. 

The unperturbed separatrix is given by Aq~ = T 0 x F u . It is foliated by the one-parameter 
family of integral curves 

r i 0 : {%{t),Po(t),t -t 0 ), (2.14) 

to £ F u being the parameter. Let 7 : [0,1] —> C be a closed path in the complex plane with 
7 ( 0 ) = 7 ( 1 ) elcC. The following function on F u 

d(t 0 ) := f {H 0 , Hi}(q 0 (t),p 0 (t),t — t 0 )dt (2.15) 

J'y 

is usually called Poincare-Arnold-Melnikov integral. Here {,} is the Poisson bracket. Then 
the following result is valid: 

Theorem 4. (Ziglin) If the function d(to) has a simple zero, then for sufficiently small \ e\ 

0, the separatrix Af has a transversal self-intersection and the system \2.12 1) has no additional 
holomorphic first integral. 

It appears that there is a relation between Theorem 2 and Theorem 4. Morales-Ruiz [ 15] 
proved that, under certain assumptions, the Ziglin’s condition about the Poincare-Arnold- 
Melnikov integral can be interpreted by the fact that the Galois group of the perturbed 
variational equation along the integral curve To is non-abelian. In other words, if Poincare- 
Arnold-Melnikov integral d(to) is not identically zero, the Galois group of the perturbed 
variational equation is not abelian and the system is not integrable by means of meromorphic 
first integrals. 


3 Proof of Theorem 1 

In what follows we assume that t , qo(t ), qj(t) are complex quantities, but we keep the param¬ 
eters real. The proof goes in the following lines. For the first two cases we find particular 
solutions. Then we study the variational equation (VE) along these solutions. The first case 
is the simplest, that is why we start with it. The variational equation (VE) is reduced to a 
particular case of double confluent Heun equation, which Galois group is more or less known. 

The second case needs more steps. The identity component of the Galois group of (VE) 
is not commutative except for some discrete values of $bf- By studying higher variational 
equations we find a logarithmic term in solutions of (VE 2 ) and (VE 3 ) when g^F 7 ^ 0, which 
implies non commutativity of the identity component of Gal(L 2 /K) (Gal(L 3 /K)) and hence, 
non-integrability of our Hamiltonian system. 

For the third case we use a perturbational technique which is still related to the Differential 
Galois approach. We study the Poincare-Arnold-Melnikov integral in order to show that a 
complex separatrix self-intersects. 
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3.1 The case Cq = 0, Cj ^ 0. 

In this case the Hamiltonian (1 1. 8 j) becomes 

2 -i N t N f 

Po , 1 


N f 


N f 


H = ^ + ^J2p‘j + w o9o + J] - g BF q 2 0 J2 q 2 j - j + 

i b 


(3,1) 


1 1 1 
The equations corresponding to the Hamiltonian (13.ip are 


N f 


Qo = Po, Po = -2cu 0 go + 2go + ^9bfQo ^ bA 

i 

C 2 

Qj = Pi, Pi = - 2uj iQi + 2^BFgo?i + i = 1, ■ ■ ■, Nf. 

Pi 

Proposition 1. The system \S. 2\) has a particular solution of the form 

qo=Po = 0 , 

sinh(2iy2a^), Pj = q j: j = l,...,N f . 


(3.2) 


2 C) 

Pi = 


\Z 2uj i 


(3.3) 


Proof. We put qo = po = 0 in ([372]) . The general solution of the system with respect to 
( Qi,Pi),3 = I,---, N f is 


o - 2 = A. + 

J 2 Wj 


C 2 h 2 . 

2ot- ~ AuP Sinh ~ to), Pi = 4i, J = 1, • • •, N f , (3.4) 

J 3 


here hj are arbitrary constants. Then we set hj = 0 and to = 0 to obtain our particular 
solution. 

□ 

Denote the variations by £o = dq 0 and go = dp 0 . It is easy to be seen that the (NVE) are 
written in variables £ 0 , Vo, namely 


£o = Vo, Vo = 


We rewrite (13. 5 p as a second order equation 


N f 


-2u 0 + 2g BF Qj“ 


£o- 


(3.5) 


£o + 


Nf 


c 


2u 0 - 2g BF V' J— sinh( 2 iy / 2 o)~f) 
i V 2uj i 


£,o = 0 . 


(3.6) 


The study of the identity component of the Galois group of (j3.6[) is a difficult task. That is 

2 _ 

why we assume that all ujj are equal. We put Uj = j = 1,... ,Nf. Then we get a variant 
of Mathieu equation 

lo T [-'ll + sinh(2i(nf)] £o = 0, (3-7) 


where 


Nf 


Zj 

A\ = 2cn 0 , Bi = - g BF > Cj. 

uj 


(3.8) 
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Since Cj are constants of integration, we can always assume that Cj 7 ^ 0. 

Next, by changing the independent variable x = e 2luJt we get an algebraic version of (13.71) 


£>' + -to + 


X 


BAB 

X X 2 X 3 


£o = 0, 


(3.9) 


where ' = -C, A = B = — It is obvious that when B — 0 this equation becomes an 

Euler equation which is solvable. Further, we reduce (13.91) to the standard form by putting 


V = V^o, 

y" = r{x)y, r(x) = (3.10) 

The equation (13. 10p is a particular case of double confluent Heun equation. For this equation 
the points 0 and oo are irregular singular ones and one natural way to study the Galois group 
is the Kovacic algorithm. This is done by A. Duval and M. Loday-Richaud in [ 6 ] p.237. We 
just apply their result which simply says that if B ^ 0 the Galois group of (I3.10p is SL(2, C). 
In our case 


B 


A If 

9 bf yb 
4 ca 3 ^ h 


which means that under the assumption Cj ^ 0 


B — 0 <S=7 < 7 bf — 0, 


that is, the identity component of the Galois group is noncommutative if #bf 7 ^ 0. Therefore, 
by Theorem 1 the Hamiltonian system (13. ip is non-integrable unless #bf = 0. This finishes 
the proof of this part of Theorem 1. 

Remark 2. Let us note that in EM a systematic procedure is presented, called 
Hamiltonian Algebrization, which transforms second order linear differential equations with 
non-rational coefficients into differential equations with rational coefficients. As an example, 
the Mathieu equation is considered, see for instance, section 2.1 in pjj. The conclusion is the 
same: the Mathieu equation is not integrable for B ^ 0. 


3.2 The case Co ^ 0, Cj = 0. 

Let us find a particular solution first. 

Proposition 2. The Hamiltonian system generated by the Hamiltonian H1.8\) with Cj = 0 
has a particular solution in the form 

2 - 

Qo{t) = 3^0 + p(t; 02,03), Po(t) = qo(t), qj=pj = 0, j = l (3.11) 

where pit: g 2 , g 3 ) is the Weierstrass elliptic function satisfying 

T : v 2 = 4u 3 - g 2 v - g 3 (3.12) 

with g 2 = yWg — 4 h, g 3 = 4 Cq — |u;o h + and h is level of the Hamiltonian \1.8\) . chosen 

so that A — g\ — 27 g 2 ^ 0. 














(3.13) 


Proof. We put qj — pj — 0, j — 1,..., Nf (recall Cj = 0) in (II. 8ft to obtain 


H _ P 2 o , .. n 2 Qo , Cl _ h 
ff- T + ^9b- T +2^-2- 


We rewrite this expression in the form 


C 2 

Qo = — 2o;o + Qo - T + h- 

% 

Then denoting u — qf } and also u = v + |o; 0 we obtain the general solution of (j3.14j) 


(3.14) 


o Zi 

%(t) = + p(t - «o;92,93), Po(t) = ®(t). 


(3.15) 


We set t 0 = 0 to get the desired result. 

□ 

Next we write the variational equations (VE) along the particular solution (13. lip . Denote 
£o — dq 0 , i] 0 = dp 0 , = dqj, rjj = dpj. Then the (VE) can be written as 


/ 3C 2 \ 

io = Vo, Vo = 2cuo + 6<?o (t) — _ 4 ^ J ^ 0 ) (3.16) 

ij = Vj, Vj = (-2 Uj + 2g BF qZ(t)) £j, j = l,...,N f . (3.17) 


The equation (13.1611 forms the tangent part of (VE) and the equations (13. 1 7jl form the normal 
part of (VE), actually (NVE). It is seen from (j3. 17[) that (NVE) splits into a system of Nf 
independent equations (NVEj), j = 1 Hence, (NVE) is integrable if, and only if, 

each of (NVEj) is integrable. In other words, the identity component of the Galois group of 
(NVE) is solvable (commutative) if, and only if, each of identity components of the Galois 
groups of the (NVEj) is solvable (commutative). Therefore, it is enough to study one of them. 
Let us write (NVE^) for certain particular j as a second order equation 

ij + (2 Uj - 2g BF ql(t)) ^ = 0. (3.18) 


Taking into account the particular solution (13. lip the Eq. (I3.18p is a Lame equation 


ij + ( 2c Oj — -gBFUJo — 2g B Fp(t) J ^ — 0. 


(3.19) 


It can be proven that if g B p ^ n ^ 2 +1 ^ , n G Z the monodromy group of (I3.19P is not Abelian 
(see e.g. m- Since the equation (13. 19ft is a Fuchsian one, the monodromy group generates 
the differential Galois group, and hence the Galois group is not abelian. Then due to Theorem 
[2] the Hamiltonian system with the Hamiltonian (1 1.8 p is non integrable in Liouville sense. 

Further, we study the tangential part of the (VE) - Eq. (13.161) . The theory gives that its 
Galois group is solvable. In fact, we have 


Proposition 3. The Galois group of /IV 1 61) is abelian. 

Proof. It is well known that the system (I3.16P has a particular solution (£o,i,£o,i) = 
(po(t),Po(t))■ The other solution is obtained via D’Alembert’s formula 


£ 0,2 


£o,i 



dr 

(W 5 ' 
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Denote the coefficient field of (j3.16[) by K = C (p(f), p'(t)). This field is isomorphic to the 
field of meromorphic functions «M(T) on T. 

It can be seen from the obtained solutions that one part of them lie in a quadratic extension 
of the field K and the another part is obtained with single quadrature of the elements of this 
extension. Therefore the Galois group of (13.161) acts in the following way: a G Gal(L/K), 
or(f 0 ,i) = f 0 ,i and = £ 0,2 + D)£o,i 5 a o G C. Let S(t) is the fundamental matrix of (13.161) 

~ ( ^°’ 2 ^ 

V £0,2 £0,1 / 


Then a G Gal(L/K ) can be represented by the matrix R Uo , crH(f) 


1 0 
1 


It is clear that the group 


1 0 
z/ 0 1 


is abelian. 


E(t)R U0 , where R Uo = 


□ 


So far we have shown that if gsF 7 ^ , n G Z the identity component of the Galois 

group of (NVE) is not abelian and hence the Hamiltonian system under consideration is 
non-integrable. 

Now, let us consider the case when 


9 BF 


n{n + 1) 


2 


J 


n G Z. 


Then every equation (I3.19j) is a Larne equation in Weierstrass form 


(3.20) 


~ H n + + Bj] = 0, (3.21) 

where Bj = ^LOon(n + 1) — 2c Oj. The cases for which the Lame equation (13.211) is solvable are 
well known: 

(i) The Lame and Hermite solutions. In this case n G Z and g -2 , .9.3, B are arbitrary 
parameters; 

(ii) The Brioschi-Halphen-Crowford solutions. Here m:=n + l/2GN and the parameters 
g 2 ,g 3 ,B must satisfy an algebraic equation. 

(iii) The Baldassarri solutions. Now n + 1/2 G |Z U |Z U |Z \ Z with additional algebraic 
relations between the other parameters. 

Note that in the case (i) the identity component of the Galois group G° is of the form 
^ and in the cases (ii) and (iii) G° = id (G is finite). And these are the all cases when 
the Lame equation is integrable. 

Therefore, together with the result of Proposition [2] we have that the identity component 
of Galois group of the (VE) is represented by the block-diagonal matrices of the kind 


1 0 

V 3 1 


/1 

0 

0 

0 

0 

0 \ 

D) 

1 

0 

0 " 

' 0 

0 

0 

0 

1 

0 

0 

0 

0 

0 


1 ' 

' 0 

0 

0 

0 

0 

0 

1 

0 

^0 

0 

0 

0 " 

VNf 
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and it is clearly commutative. 

The integrability of Hamiltonian systems with two degrees of freedom which (NVE) are 
Lame equations is studied in [TTJ ,12] • We summarize the facts and the result (Theorem 5), 
that gives necessary conditions for integrability in the Appendix. Since in our case the (NVE) 
splits into a system of Nf equations, the result for two degrees of freedom can be applied. 

The potential <p(qo) is obtained from fj3.14|) . Denote a(t,h ) := n(n + 1 )p(t) + Bj. We 
calculate the coefficients of the polynomial P(a, h ) (compare with (A.4) and (IA.5D ) 

P(a, /i) — (cq T /u 2 2 )ct 3 T (b\ T h& 2 )o£ T (c^ T /ic 2 )ct T (<A V hd 2 ). 

In our case these coefficients are: 

cii = 7 4 ,. , 02 = 0 , b\ — — T—n \ , b-2 = 0 

L n{n-\-l) t * ’ L n{n+ 1 ) ’ ^ 

i o r2 

C 1 = iipm) “ fu$n(n + 1), c 2 = 4n(n + l), (3.22) 

d\ = Ybj^n(n + 1) - B] - n 2 (n + l) 2 (4 C$ + \ fwg) , d 2 = 8 n(n + 1 )uj. 

Now we are ready to apply Theorem 5 (see the Appendix). 

The condition 3 is not fulfilled: c 2 7 ^ 0 and c 2 fci — 3a\d 2 = —32ui 0 n(n + 1), which is 
nonzero by the assumption that uJo,iOj are positive numbers, made in the very beginning. In 
particular, there are no Baldassarri solutions. 

We proceed with the cases of the condition 2. In the case 2.1, m — 1, b\ — 0 is equivalent 
to 

Bj = 0, or, ujj = cjo/4, j = (3.23) 

If for some j Bj 7 ^ 0 then the system is not integrable for this m. We will consider the case 
when all Bj = 0 in what follows. 

The case 2.2 m = 2 does not occur here since c 2 7 ^ 0. 

In the case 2.3, m — 3 the necessary conditions 

16 aid 2 + 116iC 2 = 0, 1024a 2 di + 704ai6iCi + 455 3 = 0 


yield correspondingly 

32 

Bj = —Uj, (55(n 0 = 28a;,-), 7 3 C 2 = 72a; 3 . 


(3.24) 


If any of the above conditions is violated, then the system is non-integrable for this m. We 
will consider the case when relations (13.241) are valid for all j in what follows. 

Finally, the case 2.m, m > 3 does not occur here since c 2 7 ^ 0 and d 2 7 ^ 0. 

In order to resolve the condition 1, the case 2.1 with (13.231) and the case 2.3 with (13.241) 
of the condition 2 in the Theorem 5 we need to study the Galois groups of higher variational 
equations and to apply Theorem 3. To compute higher variations we put 


qo = 

+ 

II 

efo 1} 

+ 

^ 2) 

+ e^ 3) 

+ 

Po = 

+ 

0 

II 


+ 

£ 2 Vo ] 

+ £ 3 Vo ] 

+ 

dj = 

= 0 + 

*j 1J 

+ 



+ 

Pi = 

= 0 + 

£r l? 

+ 

£ 2 r]f^ 

+ e 3 rjf^ 

+ 
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(3.25) 
















and substitute these expressions into the original Hamiltonian system. Comparing the terms 
with the same order in e, we get consecutively the variational equations up to order 3. 

The first variational equation is 


£o 1} = Vo\ Vo’ = (^-2w 0 + 6 q£ 

C (1) = vf\ V? = (-2 Uj + 2g BF ql)&\ j = 1,..., N f , 


,(!) 


•W 


■ —2 


3 C? 


9o 


0 \ M) 
so > 


*3 O’ >3 y J 

but of course we know it (see (13.161 (13.171) 1. For the second variational equation we have 


(3.26) 

(3.27) 


A = A, rtf’ = ^- 2 W „ + 6 ft - ZZj + hr, 

ij 2) = vf\ vf ] = (-2 Uj + 2g BF ql)if ] + I<f\ j = 1 


The third variational equation is 


Here 


7(3) 

SO - 

(3) 

= Vo , 

A = 

( „ 3Cq 

2ta 0 + 6 g5 4 

7(3) _ 
S j 

(3) 

= v)\ 

V? = 

-2c oj + 2(7bfQ , o)£] 3 ' ) 

Af> = 

^9bfQo y^(£^) 

2 + 6?o(£o^) 2 + 6 Cq 

Kf = 


j = 1, • • •, Nf, 

Af> = 

2 gBF 

kEW+tf'DfS 11 ) 


‘(3) 


'(3) 


(£. 


'(3) 


( 1)\2 


—ic; ? 

% 




t(l)t(2) 

.0 so 


$ [ioi ^ 1 ’) 3 - mAA 
% 1 


Jlj 3) — 2g BF 


Then, in our notation from Section 2, we have 


Kjqyp + 2* (aA + dhfll. j = 


(3.28) 

(3.29) 

(3.30) 

(3.31) 


(3.32) 


h = 
h = 


0 , Kq , 0 , Aq ,..., 0 , K 


N f 


0 K^ 0 A'® 0 K ^ 

u >- fY 0 ) • • • > u ; ly N f 


1 T 


(3.33) 


First, we have to solve (VEx). Let £ 0*1 >£02 f wo linearly independent solutions of (13.26j) 

with Wronskian equal to unity, i.e., £oq£og - £oq£og = 1- Similarly, ^, £]'9 are linearly 
independent solutions of (13. 27[) with Wronskian equal to unity. Then the fundamental matrix 
X(t) of (13.261) . (13.27ft and its inverse have the block-diagonal form 
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[til til oo o o \ 

til til oo'" oo 


m = 


0 0 
0 0 


,(D 

H 1 

t(i) 

s, i 


,(D 

Si, 2 

7« 

S?,2 


0 0 
0 0 


x-\t) 


0 0 


0 0 


yd) 

SjV/,1 


yd) 

SjV f ,2 


V 

0 0 

0 0 


t(i) 7d) 
SJV/,1 SiV/,2 

/ 


/ 7 (1) 

SO,2 

_yd) 
so, 2 

0 

0 

0 

0 

\ 

_yd) 

S0,1 

yd) 

S0,1 

0 

0 

0 

0 


0 

0 

yd) ' ' 
Sj,2 

lyd) 
sj ,2 

0 

0 


0 

0 

_7(D 

Sj,l 

td) 

H 1 

0 

0 


0 

0 

0 

0 

7(i) 

SVf ,2 

••• _yd) 

SAT^l 

Ltd) 

SiVj,, 

2 

0 

0 

0 

0 

yd) 

SlVy,! 

/ 


(3.34) 


(3.35) 


The first variational equations (VE]) (13.261) . (j3. 27[) have a singular point at t — 0 (the 
pole of p(t)). We calculate the expansion of the solutions of variational equations around the 
point t = 0. Note that 


- f .s _ 1 ,Uo + ,( 92 ul 

9o(i) “ 7 + T* + (« _ Is 


r + • • • • 


(3.36) 


Here and further dots denote the higher order terms with respect to t. In a neighborhood of 
t = 0 we have the following expansions for the solutions of the tangential part of (VEi) Eq. 
(1T26D 


(1 ) = _ ~u _ 

Wl 4-2 Q 


(U 0 


3#2 
40 


UJn 


r +... 


(!) _ ^ W 0 5 + 

4o ’ 2 “ 5 + 35* + "" 


(3.37) 


t 2 3 V 40 6 

Now, we suppose that #bf = "L +1) 7 ^ 0. 

First, let us consider the condition 1: a 3 = , 4 

’ n(n+l) ' 

(i). We take n — 1 for simplicity, but we keep writing ^bf instead 1. In the vicinity t = 0 we 


n G Z, i.e., the Lame and Hermite case 


have the following expansions for the solutions LLLL °f 113.271) (n = 1) 


^ = 7 + ^+1 + = + 




9 2 




(i) 


a j_. 

30 


(3.38) 


where Bj = 2 Uj — 4a;o/3. 

There are no logarithms in the expansions around t = 0 of the local solutions of the second 
variational equation (VE 2 ). 

We will show that a logarithmic term appears in a local solution of (VE 3 ). For this 
purpose, it is enough to show that at least one component of W _1 / 3 has a nonzero residue at 
t — 0, see formulae (I2.7p . (I2.8|) . We calculate j-th component j — 1,..., Nf of X f _1 / 3 , which 
looks like 




(3.39) 
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We take 


t(i) _ t(i) c(i) _ t(i) 

SO — SO,2) Sj ~ Sj,l- 


With this choice we find 


J2) _ gBpNf 1 _ Wo 

* 0,1 - t 2+ 2 t 3 + 


^( 2 ) 9BFNf 1 : 

4°, 2 - 2 t + • • • 


and 


y( 2 ) — - 4- ^-t 4- yw — _ 4 . 

Sj,l - , + o 1 + ' ' ' ’ Sj,2 - o +- 


1 Bj 


( 2 ) 


Taking the first term in (13.391) . namely 


M3 = -$*?> = -4‘»2 Jbf [«««« + 2* (^»if+of) 


(3.40) 


(3.41) 


(3.42) 


/n\ /o\ /o\ /o\ 

with the choice (I3.40p and ^ = Q ,2 and £• ' = we can see that /i has a simple pole at 
t = 0 with residue —2g FF Nf/3, which is non-zero. Therefore, the identity component of the 
Galois group of (VE 3 ) is not commutative and hence, in this case, the Hamiltonian system 
QL 8 P is not integrable due to Theorem 3. 

Similarly, for n — 2 and #bf = 3 we have the following expansions for the solutions ^ 

offlMID 


® = h~T + °W 


yd) ^ i B-jt'' 

5 + 70 


+ ..., 


(3.43) 


where Bj = 4co 0 ~ 2uij(n = 2). Let us study first the expansions of the local solutions of the 
(VE 2 ) around t — 0. The calculation of /i 2 = K^ 2> with ^q 1 - 1 = £o*i,£^ = gi yes 


^ 2 


12 45, | B]- fg 2 B,(g 2 - l B j) 

C i? t 3 + t {> 


Since g 2 depends on h, which is arbitrary provided A = g 3 — 27 g\ 7^ 0, the only possibility for 
the residue of /z 2 to be zero is Bj = 0 or ojj = 2cu 0 . If at least exists one uj v such that Bj 7 ^ 0 , 
then a logarithm appears in the solutions of (VE 2 ) around t = 0. 

We proceed with the case when all Bj = 0, or equivalently, c Oj = 2u 0 ,j = 1 
Choosing 


t(i) _ t(i) _ 1 _ ^0 , yd) _ yt 1 ) _ f_ 

so S0,1 j2 3 ’ S j Sj,2 g 


(3.44) 


we find that 


and 


y(2) 1,1 ^0 ^0 , 

^°d “ t 3 + t 2 5 t 3 + 


y(2) _ 1 <^0 

^°’ 2 ~~ t 3 5 t 






Taking the second term in (13.391) /i 3 = with the choice (13.441) and £0 21 = ^ 2 , ^ = 

one can see that /J 3 has a simple pole with a residue —cuo72/25, which is nonzero since oj 0 7 ^ 0 
by assumption. 

In either of the cases above, the identity component of the Galois group of (VE 2 ) or (VE 3 ) 
is not commutative and the Hamiltonian system (11.81) is not integrable due to Theorem 3. 
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Next we consider the case 2.1 with (13.231) . Here n = \ and #bf = §• We take 

*(i) _ t(i) t(i) _ A i) 

SO >0,1 ’ S7,2 ’ 

There are no logarithms in the expansions around 1 = 0 of the local solutions of the second 
variational equation (VE 2 ) due to (13.2311 . With the above choice, we find that 

«S 2) = eg = p + ow, = e' 2 ,’ = 2= - ^+o(« 3 ' 2 ). 

Then the first term in (13.39j) has the following expansion around t — 0 


/A 


A 1 ) 7V'( 3 ) 

'S',2 j 


3 2 
8 l 2 


~3t~ + 


that is, n 3 has a pole at t — 0 with non-zero residue Therefore, the identity component 
of the Galois group of (VE3) is not abelian and hence, in this case, the Hamiltonian system 
(11. 8p is not integrable due to Theorem 3. 

Finally, we consider the case 2.3 with (13.2411 . Here n — § and #bf = y- We take 


*(1) _ Ai) 

SO — S0,1) 


t(!) _ t(l) 

S S, 1 ■ 


There are no logarithms in the expansions around 1 = 0 of the local solutions of the second 
variational equation (VE 2 ) due to (13.24ft . With the above choice, we find that 


_ A2) _ 51V/ 1 _ iV/^o 

?0 ’ 2 1441 4 1 3 2241 2 


(Uo 

51 


+ o(0, 


ef = = * 7/2 (i+o<* 2 )+ 1 - 7/2 (4 - i (2 + ■ ■ ■) ■ 

Again the first term in (I3.39P has the expansion around 1 = 0 


_ A 1 ) zy( 3 ) _ 

^3 - t 3 ,2 K 3 ~ 


175 N f 35 5Nf0j 0 7u 0 0 

17281 4 31 s 5761 2 121 1 


that is, /i,3 has a non-zero residue at 1 = 0. Therefore, the identity component of the 
Galois group of (VE 3 ) is not abelian and hence, in this case, the Hamiltonian system (1 1.8 p is 
not integrable due to Theorem 3. 

Remark 3. For arbitrary neZ in $bf = " ( ' i 2 +1 ' > one needs to know the exact coefficients 
in expansions of the Lame solutions of (13.27ft and eventually the expansions of the higher vari¬ 
ations. The formulas are quite involved. However, it is unlikely that the system is integrable 
for some n > 2. 

This finishes the proof of this part of Theorem 1. 


3.3 The case Cq ^ 0, C\ ^ 0. 


Here we consider the Hamiltonian (11. 8ft only for two degrees of freedom (see comments in the 
next section) 

77 = tt + UoQo ~ 7T + tS + 7T + w i0i + tS - 3bf9o9i- (3-45) 


2 q { 


0 


2 q( 
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(3.46) 


Denote e := pbf and assume that e is small enough. We can rewrite (I3.45P as 

H = H 0 + eH 1 , 

where 


n 2 n 4 C 2 n 2 C 12 

TT - Po _L , , n 2 I , Pi , , , 2 I L 1 

no — — + ^0% -— + 77^- + TT + 0J\q t + —2, 


The unperturbed system (e = 0) is separable. 


2 q( 


Hi = -q 2 0 ql 


(3.47) 


3 , c 2 


Qo — Poi Po — — 2 u 0 q 0 + 2 q 0 + 

c 2 

qi=Pi, P i =-2wigi +-4 

«1 


0 

«iT’ 


(3.48) 

(3.49) 


From the proof of Proposition 2 the general solution of (13.481) is found in (13.151) . From the 
proof of Proposition 1 the general solution of (13.49(1 is 


2 _ 

qi -%7, + 


cl hi 


2u>i 4uf 


sinh 2iy/2ui(t — to), pi = q\. 


(3.50) 


First, we put the Hamiltonian H in the context of the theory recalled in Section 2. It is 
assumed that at this point the variables are real. We introduce action-angle variables (I,p), 
so that H 0 = H 0 (q 0 ,p 0 , /). To do so, we need to find a generating function S(I, q\) : 

dS 


such that 


S(I,qi) dS 

ipi,qi) ^ Pi = Q t -> 


TP + upql + = hi -A hi(J) := /. 


2 <?r 


(3.51) 


Note that the real ovals for the curve (pi, gi) in (13.5ip exist for h\ > C] . Then the formula 

V C 1 

(I3.50P becomes 






The generating function S' can be found explicitly, but we do not need it, we just set 

, _ Pi , 2 . _ fdqi 


(3.52) 


(3.53) 


Note that dl A dp = dpiAdqi, p is multivalued, but p — 1, that is, t and p are interchangeable. 
Next, we fix / to an arbitrary constant greater than Cl and again consider t, qo(t),p 0 {t) 

V C 1 

as complex variables. Our system becomes an one-and-a-half degrees of freedom system with 
a Hamiltonian H — Hq + eH\, where 


tt _ Po , 2 <?o , Co j jj _ 2(1 

~ Y +UJ ° q o~Y + 2q 2+1, Hl ~ ~ q ° 


1 2 


10 


Cl 

2uji \l 4cuf 2u>i 


sin 2\/2u/i(t — t 0 ) . (3.54) 
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We need to find a separatrix in the dynamics of ( qo,Po )• Denote h = h — I and g- 2 = 
yC^Q—4/z, g 3 = 4Cq — Iwoh+yWo (compare with the corresponding formulas in the Proposition 
m ■ Let h* be the biggest real root of 


A (h) = g\ — 27gl = —64 ^ h 3 — co^h 2 — QC^uJoh + 8Co<^o + ^j-Co^ = 0- (3.55) 

Assume that 4u j$ — 3 h* > 0. Further, we denote 

a := 3h • > Q 

Then the unperturbed system (13.541) has a separatrix 

2 3 

r o : q%(t) = -u 0 + a+ Q Po(t) = 9o(*)- (3-56) 

o smh (V3ar) 

The perturbed variational equation (PVE) of (j3.54[) along r 4o is given by (see [71115]) 


d (^\ ( 7/()po/) 0 #i, P o \ /s \ 

^ I V I = I —Ho,qoqo —Ho,qopo ~Hi qo I I Tj I , 


(3.57) 


where all coefficients are restricted to r to . In order to study the Galois group of (PVE) we 
fix the coefficient held K in (I3.57p . From the expressions for the separatrix (I3.56P and the 
perturbation H 3 (I3.54p 

K := c( e ^“ t ,e 2V5 ^ it ). 

Then, to obtain the fundamental matrix of (13.57j) a quadrature is needed, namely 5 — 8(t) — 
f H ^ 0P0 dt (see [15J for details). In our case S — f V equals 

0,P0 | r p 0>. ) 


5 = — —y(— ° sinh(\/3at) cosh 3 (\/3 at) H— ° )L_ sinh(\/3at) cosh(\/3at) 


(3a ) 3 12x730 


8\/3a 
2 c 0 q -I - 12 a 


3\/3 a 


tanh(x/3af) + 


26a;o + 99a 


t). 


It is clear that 8 — 8(t) is uniform and 8 ^ K. Then, the Picard-Vessiot extension of (I3.57p 
is Li = K{8) = C(e x/3 “ f , e 2 '^ 2 ^ lt , t). It remains to hnd d(t 0 ). Let 7 be a loop around the pole 
t — 0. Then simple calculations give that the Poincare-Arnold-Melnikov integral is 


d(t 0 ) 


{H 0 ,Hi}(q 0 (t),p 0 (t),t 


t 0 )dt = Y2ma\/2(jJi 


1 2 


C 2 , - 

-sin 2 \/ 2 coito- 

2u 1 


(3.58) 


It is seen that d(t 0 ) has simple zeroes and by Theorem 4, the perturbed separatrix self¬ 
intersects transversally. Also since d(t 0 ) is not identically zero, the Galois group of the 
perturbed variational equation is not abelian ng. Hence, when e = ^bf 7 ^ 0 sufficiently 
small, there is no additional meromorhic Erst integral. This finishes the proof of this part and 
therefore, the proof of the Theorem 1. 
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4 Concluding Remarks 

In this paper we use variational equations to obtain a necessary and sufficient condition for 
integrability of a system which describes the stationary solutions in the time dependent mean 
held equations of Bose-Fermi mixture. Here we make some remarks. 

We start with some restrictions to our methods. In subsection 3.1 we don’t know how to 
study the Galois group of a second order linear equation with quasi-periodic coefficient, that 
is why we assume that all ujj are equal. It is an open problem to develop a Picard-Vessiot 
Theory for the coefficient held K = C(e aiX , ..., e arnX ) with an, ... , a m £ C and to relate this 
result with the integrability of the corresponding linear equation, see [18], p. 408. 

In 3.2 we consider n = 1 and n = 2 only by technical reasons. It is highly unlikely that the 
system is integrable for n > 2, n £ Z, which is justihed by the result in 3.3. We notice that 
the non-integrability result obtained in this case are also valid for the limiting case Co = 0 
and Cj — 0, j = 1,..., Nf. 

In the general case Co ^ 0 and Cj ^ 0, j = 1 ,Nf (VE) does not split in nice way 
as in the previous cases. Because of this reason, we consider the system (11.81) with two 
degrees of freedom. Even then, studying the Galois group of (NVE) is not so simple due 
to a number of parameters. That is why we use a perturbational approach, which is still 
related to the Differential Galois approach. Furthermore, this approach gives a dynamical 
meaning to the algebraic obstructions to integrability. Note that, the using Poincare-Arnold- 
Melnikov integrals in more degrees of freedom for real Hamiltonian systems needs certain 
KAM-conditions. 

The above results allow us to think that the system (11.81) is not integrable unless #bf = 0. 
Moreover, the formulas (13.41) and (13.151) give the general solution to the separable system 
(S'bf = 0). 

Acknowledgements O. C. acknowledges funding from Bulgarian NSF Grant DDVU 
02/90. 
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A Necessary conditions for integrability of Hamiltonian 
systems which have (NVE) of Lame type 


In this appendix we recall some facts concerning the integrability of Hamiltonian systems 
with two degrees of freedom, an invariant plane and which (NVE) are of Lame type. More 
details can be found in mm • In our case the (NVE) splits into a system of Nf equations 
of Lame type, and therefore, these arguments can be applied. 

Classically the Lame equation is written in the form 

f- [n(n+ l)p(t) A = 0, (A.l) 

where p(t) is the Weierstrass function with invariants g 2 and 93 , satisfying v 2 = 4u 3 — g 2 v — . 9.3 
with A = g\ — 27 . 9.3 7 ^ 0. 

The known (mutually exclusive) cases of closed form solutions of (IA.1D are: 

(i) The Lame and Hermite solutions. In this case n E Z and g 2 ,g 3 ,B are arbitrary 
parameters; 

(ii) The Brioschi-Halphen-Crowford solutions. Here m:=n + l/2eN and the parameters 
g 2: g 3 ,B must satisfy an algebraic equation. 

(iii) The Baldassarri solutions. Now n + 1/2 6 |Z U |Z U |Z \ Z with additional algebraic 
relations between the other parameters. 

Note that in the case (i) the identity component of the Galois group G° is of the form 

^ ^ and in the cases (ii) and (iii) G° = id (G is finite). And these are the all cases when 

the Lame equation is integrable. 

Now consider a natural two degrees of freedom Hamiltonian 

H = - (jp\ A P2) A V(qi, 92), (A. 2 ) 

qj(t) E C,pj(t) = qj,j = 1, 2. We assume that there exists a family of solutions of the form 
La 92 =7*2 = 0 , q 1 = q 1 (t,h), pi(t,h) = qi(t,h) 


and gi (t, h ) is a solution of 


The (NVE) along T h is 


1 

2 


9i + ^(9i) = h, 


h G M. 


f = 0, 


(A.3) 


where a(t, h) = a(qi(t, h)) is such that (IA.3[) is of type (lA.Ilh 

In [Hi I2J the type of the potentials V with this property are obtained as well as the 
necessary conditions for the integrability of the Hamiltonian systems with the Hamiltonian 
(I A. 2 D . In order to formulate the result we need certain additional quantities. 

Since a(t, h) depends linearly on p(t), then a 2 is a cubic polynomial in a, depending also 
in h, namely 

a 2 := P(a, h ) = P\(a) + hP 2 (a). (A.4) 


The following coefficients are introduced 


P{oi , h) — ( cl\ T ha 2 )a 3 ffi (bi T hb^ct 2 T (c^ A hc 2 ^Oi A (d\ A hd 2 ). (A.5) 


Now we are ready to give the corresponding result. Note that the following Theorem gives 
necessary conditions only from the analysis of the first variational equation. 
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Theorem 5. (Theorem 6.2 Assume that a natural Hamiltonian system has (NVE) of 

Lame type, associated to the family of solutions lying on the plane q 2 = 0 and parametrized 
by the energy h. Then, a necessary conditions for integrability is that the related polynomials 
Pi and P 2 satisfy a 2 = 0, and one of the following conditions holds: 

1- ai = for some n G N; 

2. a\ = -j^r—j- for some m G N. Then, assuming the conjecture above is true, one should have 

b 2 = 0 and we should be in one of the following cases: 

2.1) m — 1 and b\ = 0 , 

2.2) m = 2 and c 2 = 0 , 16aiCi + 3 b\ = 0 , 

2.3) m = 3 and 16a±d 2 + llbic 2 = 0, 1024a^(ii + 704ai&iCi + A5bf = 0, 

2.m) m > 3. Then, we should have b\ — 0 and, furthermore, either c\ = c 2 = 0 if m is 

congruent with 1, 2,4 or 5 modulo 6, or d\ = d 2 = 0 if m is odd; 

3. a i = n (n+ i) wn + 1/2 G |Z U -JL U |Z \ Z, b 2 = 0 and either c 2 = 0, b\ — 3aiCi = 0 or 

C 2&1 — 3 aid 2 = 0 , 26‘J — 9ai&iCi + 27a\di = 0 . 

It is clear that the condition 1. in the above Theorem gives the Lame and Hermite solutions 

(i), the condition 2.- the Brioschi-Halphen-Crowford solutions (ii), and the condition 3. - the 
Baldassarri solutions (iii). 
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